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Abstract 

Analytical  solutions  play  a  very  important  role  in  the  validation  of  numerical  codes.  However,  exact  analytical 
solutions  involving  optimal  design  of  transiently  loaded  multilayered  structures,  are  rare  in  the  literature. 

In  this  paper,  we  solve  an  optimal  design  problem  involving  wave  propagation  in  a  two-layered  elastic  strip  subjected 
to  transient  loading.  We  obtain  explicit  formulas  for  the  stress  in  each  layer  using  the  method  of  characteristics,  and 
then  use  these  results  to  identify  the  designs  that  provide  the  smallest  stress  amplitude.  The  derived  analytical  results  are 
then  successfully  used  to  validate  a  previously  developed,  hybrid  computational  optimization  software. 

©  2003  Elsevier  Ltd.  All  rights  reserved. 

Keywords:  Wave  propagation;  Goupillaud-type  layered  media;  Structural  optimization;  Analytical  solutions;  Validation  of  numerical 
codes 


1.  Prior  work 

The  study  of  wave  propagation  in  layered  media  is  important  in  a  number  of  diverse  disciplines  such  as 
seismology,  electromagnetics,  optics,  acoustics,  et  cetera.  An  exhaustive  review  of  the  literature  on  the 
subject  is  beyond  the  scope  of  this  paper,  and  the  reader  is  referred  to  the  classic  treatise  by  Brekhovskikh 
(1960)  for  a  review  of  the  mathematical  framework  and  physical  phenomena  related  to  wave  propagation  in 
elastic  and  electromagnetic  media.  The  different  approaches  to  the  design  of  layered  media  for  optimal 
reflection/transmission  characteristics,  can  be  broadly  classified  into  either  frequency  domain  or  time  do¬ 
main  methods.  Some  relevant  work  in  this  area  is  included  below.  Using  frequency  domain  methods,  Hager 
et  al.  (2000)  derive  relations  needed  to  minimize  or  maximize  reflection  in  elastic  coatings  subjected  to 
steady-state,  band-limited,  time-harmonic  acoustic  waves.  Hager  and  Rostamian  (1987)  also  investigate 
conditions  necessary  to  minimize  acoustic  reflections  from  a  wall  with  a  viscoelastic  coating.  Wesolowski 
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(1995)  studies  a  problem  for  a  two-layer  interface  sandwiched  between  homogeneous  elastic  layers  sub¬ 
jected  to  time-harmonic  waves.  Gusev  (2001,  2002)  provides  extremal  relations  for  the  optimal  design  of 
homogeneous  and  inhomogeneous  layered  structures  subjected  to  time-harmonic  acoustic  waves.  Abgaryan 
and  Lyubashevskii  (1984)  develop  relations  for  active  suppression  of  time-harmonic  waves  in  a  structure 
containing  a  nondirectional  sensor  and  radiator.  Krasil’nikov  (1983)  uses  the  method  of  “slowly  varying 
amplitudes”  to  find  optimal  designs  for  a  multilayered  interference  absorber.  Bao  and  Bonnetier  (2001) 
consider  an  optimal  design  problem  of  a  periodic  structure  subjected  to  time-harmonic,  transverse  mag¬ 
netically  polarized  waves  and  establish  the  existence  of  the  optimal  designs  for  the  class  of  problems 
investigated. 

In  contrast  to  the  body  of  literature  on  time-harmonic  wave  propagation,  there  is  much  less  published 
work  on  the  optimization  of  layered  media  subjected  to  transient  loading.  Nygren  et  al.  (1999)  study  the 
problem  of  maximizing  the  efficiency  of  elastic  energy  transmission  through  a  layer  of  elastic  junctions  for 
application  to  percussive  drilling.  Konstanty  and  Santosa  (1995)  pose  a  transient  optimal  design  problem 
for  minimally  reflective  coatings  in  the  time-domain,  and  solve  the  problem  numerically  using  a  BFGS 
updated  secant  method;  their  analysis  is  limited  to  problems  with  small  interlayer  impedance  contrast.  A 
number  of  other  authors  consider  problems  related  to  the  propagation  of  stress  waves  in  discretely  layered 
or  inhomogeneous  media;  however,  these  authors  do  not  address  the  problem  of  optimization.  For  in¬ 
stance,  Ali  (1999)  uses  transmission  line  theory  to  obtain  recurrence  relation  solutions  for  acoustic  waves  in 
a  layered  medium.  Tenenbaum  and  Zindeluk  (1992)  develop  exact  algebraic  expressions  for  the  reflected 
wave  in  a  semi-infinite  layered  homogeneous  or  inhomogeneous  medium.  Lee  et  al.  (1975),  provide  some 
general  expressions  for  the  stress  jump  in  both  functionally  graded  and  discretely  layered  media,  and 
provide  error  estimates  when  a  medium  with  continuous  property  variation  is  replaced  by  a  medium 
consisting  of  a  series  of  discrete  homogeneous  layers.  Chiu  and  Erdogan  (1999),  and  Scheidler  and  Gazonas 
(2002),  solve  several  transient  wave-propagation  boundary-value  problems  for  free-fixed  and  free-free 
boundary  conditions  in  power-law  and  quadratic  inhomogeneous  one-dimensional  media,  using  Laplace 
transform  methods. 

Anfinsen  (1967)  treats  a  boundary-value  problem  similar  to  our  own,  and  finds  optimal  material 
properties  which  maximize/minimize  the  amplitude  of  the  first  transmitted  stress  wave  along  an  elastic  strip. 
Anfinsen  develops  difference  equations  and  solves  them  using  z-transform  methods.  In  addition,  Anfinsen 
determines  asymptotic  similarity  relationships  applicable  for  optimal  two-layered  designs.  We  corroborate 
Anfinsen’s  findings  which  show  that  increasing  the  impedance  ratio  of  the  layers  tends  to  minimize  the 
maximum  stress  on  a  design.  Interestingly,  despite  the  work  of  Anfinsen,  our  approach  and  results  for  the 
stress  and  optimal  design  of  a  two-layered  elastic  strip,  appear  to  be  entirely  new.  Furthermore,  we  use  our 
analytical  results  to  validate  the  hybrid  DYNA3D/GLO  optimization  software  (Gazonas  and  Randers- 
Pehrson,  2001),  previously  used  to  find  optimal  designs  to  a  class  of  impact  and  penetration  problems  using 
constitutive  behaviors  representative  of  lightweight  multilayered  armors.  Finally,  we  mention  that  Bruck 
(2000)  develops  a  time  history  profile  for  the  stress  waves  and  discovers  a  time  delay  benefit  when  using 
functionally  graded  materials.  In  our  work,  we  discover  a  similar  benefit  when  using  two-layered  optimal 
designs  other  than  the  homogeneous.  Such  designs  provide  a  time  delay  benefit  when  the  peak  values  of 
stress  are  reached. 


2.  Analytical  solutions  for  stress  wave  propagation  and  optimal  design 

We  consider  one-dimensional  wave  propagation  in  an  isotropic  elastic  strip.  The  strip  is  assumed  to  be  of 
finite  length  L,  made  of  two  layers  of  arbitrary  lengths  L\  and  L2,  so  that  L\  +  L2  =  L.  The  layer  interface  is 
located  at  the  fixed  position  x  =  Lu  where  0  <L]  <  L.  The  left  face  of  the  strip  is  subjected  to  a  stress 
loading  p,  while  the  right  face  is  fixed  (Fig.  1). 
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The  density  and  elastic  modulus  along  the  strip  are  denoted  by  the  piece-wise  constant  functions  p(x) 
and  E(x),  taking  values  p, ,  p2,  and  E\,  E2,  in  each  layer,  respectively.  In  this  work,  we  consider  longitudinal 
elastic  waves  for  a  uniaxial  strain  state  (Fig.  1),  however  our  results  are  also  valid  for  the  case  of  uniaxial 
stress.  Using  the  definitions  of  the  wave  speed  c  =  (E/p)1/2,  and  characteristic  impedance  z  =  pc,  we  relate 
to  each  layer  the  wave  speed  and  characteristic  impedance  c\,  z,  and  c2  and  z2,  respectively.  The  uniaxial 
strain  elastic  modulus  E,  is  related  to  the  Young’s  modulus  E,  through  the  expression 
E  =  E{  1  -  v) /[( 1  +  v)(l  —  2v)],  where  v  is  Poisson’s  ratio  (Meyers,  1994). 

Our  design  parameter  a  is  chosen  to  represent  the  impedance  ratio  between  layer  1  and  layer  2: 
a  =  ^  >  0.  The  transit  time  t  through  the  strip  is  given  as 

U  U 

T  = - 1 - . 

Cl  c2 


Throughout  this  paper  we  assume  a  Goupillaud-type  layered  medium,  which  is  a  medium  that  ensures  the 
same  wave  travel  time  through  each  layer,  Claerbout  (1976), 


L\  _  L2  _  t 
ci  c2  2 


(1) 


In  order  to  preserve  the  same  wave  travel  time  for  both  layers  (ci/c2  =  L\/L2  =  p),  the  layer  properties  have 
to  relate  as:  a  =  (ppl)/p2  =  El/(pE2).  In  the  special  case  of  two  layers  of  equal  length,  p  =  1  and  the 
density  and  elastic  modulus  become  proportional. 

This  problem  can  be  easily  converted  to  the  case  of  two  layers  of  equal  length  and  wave  speed  of  unity, 
by  replacing  the  spatial  variable  x  with  the  new  variable  £  =  ^  and  using  condition  ( 1 ).  Here,  c  =  c(s)  is 

the  piece-wise  constant  wave  speed  function,  taking  values  C|  and  c2  in  each  layer,  respectively.  As  a  result, 
the  wave  equation  becomes: 
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e!«_s(zS) 

Z  Qfi  0£  ’ 

where  z  =  z(£)  represents  the  characteristic  impedance  given  by  the  piece-wise  constant  function: 

fzi,  where  0<£<§, 
z^'  |  z2,  where  ^  ^  <  t. 


(2) 

(3) 


The  wave  speed  now  becomes  the  same  (unity)  in  each  layer  (c  =  ci  =  c2  =  1),  while  the  design  parameter 
a  =  3-  >  0,  as  stated  before. 

Let  t)  represent  the  value  of  the  stress  at  position  c,  time  t,  and  design  parameter  a.  Our  goal  is  to 
find  a  design  that  provides  the  smallest  stress  amplitude  during  the  wave  propagation  along  the  strip.  We 
formulate  our  optimal  design  problem  as 

(P)  inf  sup 

a>0  (K«+oo 


subject  to  the  initial/boundary-value  problem, 
tw  =  §i It  when  £  ^  , 

<  a(0,t) =  z|(0,t)  =pH(t),  u{x,t)  =  0,  (4) 

l  «(£,<>)  =  1(5,0)  =  0. 


Here,  u(x.  t)  represents  the  displacement  at  (x.  t ),  and  H(t)  represents  the  Heaviside  function.  As  demon¬ 
strated  later  in  the  paper,  solving  problem  (P),  becomes  equivalent  to  solving  the  physical  problem  with 
layers  of  unequal  lengths,  because  the  only  essential  condition  that  influences  the  stress  wave  propagation 
along  the  strip,  is  the  equal  travel  time  through  each  layer.  Therefore,  from  now  on,  any  conclusions  made 
for  the  case  of  two  layers  of  equal  length  (Fig.  2),  will  apply  to  the  general  physical  case  of  layers  of  unequal 
lengths  (Fig.  1). 

For  the  two-layered  strip  shown  in  Fig.  2,  the  propagation  of  the  stress  wave  from  layers  1  to  2,  can  be 
expressed  by  the  following  relations, 


(5) 


Here,  a  and  er+  represent  the  stress  values  behind  and  ahead  of  the  discontinuity,  while  [a]  represents  the 
stress  jump.  The  subscripts  refer  to  1  (incident),  T  (transmitted),  and  R  (reflected)  wave.  The  above  rela¬ 
tions,  see  Meyers  (1994),  involve  the  continuity  conditions  at  the  layer  interface,  and  will  be  our  main 
reference  in  calculating  the  stress  values  throughout  this  paper. 

In  the  special  case  of  a  homogeneous  design,  when  a  =  1,  combining  (5)  with  the  method  of  charac¬ 
teristics,  one  can  easily  derive  the  time  history  profile  for  the  stress.  In  the  region  ahead  and  behind  the 
stress  discontinuity  propagating  along  the  strip,  the  stress  takes  the  following  values,  respectively. 


{0  and  p ,  during  the  time  interval  (0,  r) , 

p  and  2 p,  during  the  time  intervals  {{2k—  \)x,2kx)  for  k=  1,2,.. .  ,oo, 
2 p  and  p ,  during  the  time  intervals  (2At,  {2k  +  l)r). 

This  implies  that, 

max  t)—2p, 

0<^<T,  0<f<+oo 


(6) 

(V) 


where  according  to  our  notation,  o\  (c,t)  represents  the  stress  at  position  £,  time  t,  and  design  parameter 
a  =  1.  From  here,  an  upper  bound  for  the  optimization  problem  (P)  follows, 
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Fig.  2.  Two-layered  finite  strip  under  transient  load  (transformed  coordinates). 


inf  sup  dy_(f,  t)  <  max  o\(£,t)=2p.  (8) 

a>0  0<{<t,0<«+oo  0<{<t,0<«+oo 


Based  on  the  results  given  later  in  this  section,  we  prove  that  the  reverse  inequality  in  (8)  holds,  and 
conclude  that  the  homogeneous  design  is  an  optimal  design  for  problem  (P).  We  also  notice  that  no  design 
with  design  parameter  a  <  1 ,  can  be  optimal,  since  for  such  designs  the  stress  amplitude  always  exceeds  the 
value  2p.  Indeed,  following  the  first  transmitted  wave  and  using  (5),  it  appears  that  in  the  region  behind  the 
stress  discontinuity,  the  stress  takes  values, 


£0  = 


As  the  first  transmitted  wave  reflects  at  the  fixed  boundary  located  at  £  =  t,  the  stress  doubles  and  reaches  a 
value  greater  than  2 p.  Summarizing  the  above  observations  we  have, 

sup  ^  2p  =  max  aALt),  where  0<a^l.  (9) 

0s:^T,0<r<+oo  0<{<t,0<«+Oo  v- 


In  the  rest  of  this  paper,  we  investigate  and  obtain  optimality  results  for  two-layered  designs  with  a  >  1 .  In 
Fig.  3a  and  b,  are  given  the  Lagrangian  diagrams  for  the  case  of  equal  layer  length  under  consideration,  as 
well  as  for  the  general  physical  case  of  two  layers  of  unequal  lengths  but  equal  travel  time. 

For  a  given  a  >  1,  after  applying  (5)  and  the  method  of  characteristics,  we  find  the  following  recurrence 
relation  among  the  stress  values  T^k, 


/1(°0  •  (Tx,k-i 

—  Tn,k-l)  +  Toi,k- 3) 

k  =  4,  5, ... ,  oo, 

2 

Tl-  80CP 

r  2(3a  —  \)2p 

(10) 

(a+I)A 

la’2  ~  t  ,  1  \2  ’ 
(«+  1) 

a’3  _  /  i  i  \3  ’ 

(*+  1) 

where  /1(a)  =  ^  ) .  The  recurrence  relation  given  in  (10)  is  linear,  homogeneous  and  with  constant  coeffi¬ 
cients.  Its  corresponding  characteristic  equation, 
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Fig.  3.  Lagrangian  diagram  of  stress  waves:  (a)  equal  layer  length,  (b)  unequal  layer  length. 


Tk  -  /1(a)  •  Tk~l  +  0(a)  ■  Tk~2  -  Tk~ 3  =  0, 

has  three  distinct  roots, 

{a.  -  \)  +  2y/aJ  (a  -  1)  -  2v/aJr 

y  i  =  1,  Jh  = -  , - ,  and  y3  = - — - , 

a  +  1  a  +  1 

where  J  =  yf—\.  The  general  solution  is 
T'ajr  =  A\^[  +  ^2>2  + 

with  coefficients  A\,  A2,  A 3  determined  from  the  boundary  conditions  given  in  (10), 

P  P 

A\  =  p,  A2  =  ~ ,  A  3  =  --. 

Substituting  (11)  and  (13)  into  (12),  we  obtain, 


p 

’  (a  —  1)  +  2 

k 

P 

'(a- 

1)  —  2v/oc/' 

2 

a  +  1 

2 

a  +  1 

Further  simplifications  of  (14),  lead  to  the  following  expression  for  the  term  Ty  k, 
Tx,k  =  p[  1  -  cos(fap(a))], 

where  cp(ot)  =  arctan^,  and  a  >  1.  In  fact,  the  following  relations  are  equivalent, 

2y/a 


cp  =  arctan 


a  —  1  ’ 


where  a  >  1 
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1+C0S<P  ,  A  7t 

a  =  - - ,  where  0  <cp<~.  (17) 

1  -  cos  cp  2 

From  (15),  one  can  derive  the  following  bounds  for  the  stress  amplitude  in  layer  1, 

sup  ffa(£,t)  < 2p,  where  a  5=  I .  (18) 

0<£<5,  0<r<+cx> 

These  results  indicate  that  for  all  designs  with  a  ^  1,  the  stress  amplitude  in  layer  1  never  exceeds  the  value 
2 P- 

Similarly,  using  the  continuity  conditions  at  the  layer  interface  and  the  method  of  characteristics,  we 
derive  that  the  new  stress  terms  S^k  satisfy  the  following  relation, 


Sa,k  —  Safi—  1 
=  277i  =  - 


k=  1,2 . 


or  equivalently, 

rs^  =  2(-l)*Ef=1(-l)%,  k=  1,2,.. .,00, 

1^  =  27;,!  =  ^. 

Substituting  the  expression  for  Tyk  in  (20),  we  further  obtain, 

k  k 

Sx,k  =  2(-l)  p^(-iy  -2(-l)  p^(-l)ycos(j'(jf)(a)),  (21) 

7=1  7=1 

where  as  before  q>( a)  =  arctan(2^),  and  a  >  1.  The  calculation  of  each  partial  sum  and  further  simplifi¬ 
cations  imply  that  the  term  S^k  is  given  by  the  formula, 

r  cos((2t  +  1)^)1 

cos  *4  J’  (  ’ 

Due  to  the  above  expression  for  the  terms  S^k,  the  stress  amplitude  in  layer  2  will  either  reach  or  exceed  the 
value  2 p,  for  any  given  design  with  a  >  1 .  A  more  detailed  discussion  is  provided  when  the  set  Ky  is 
introduced  in  (26).  As  a  result  we  have, 

sup  <Ta (£,t)  ^  2 p,  where  a  ^  1.  (23) 

0<f<+oo 

Combining  (18)  and  (23)  we  have  that, 

sup  a a(£,t)  ^  2 p,  where  a  ^  1.  (24) 

0  ^  ^  ^  T,  0<f<+oo 

Finally,  combining  (8),  (9)  and  (24),  we  conclude  that, 

inf  sup  =  max  =2p.  (25) 


0  ^  t,  ^  T,  0  ^  f<+oo 


O^^t,  0^f<+oo 


This  clearly  indicates  that  the  homogeneous  design  with  a  =  1,  is  an  optimal  design.  Generally,  a  design 
with  a  >  1  is  optimal  iff  the  stress  values  given  by  S%k  for  k  =  1,2,...,  oo,  do  not  exceed  2 p. 

Indeed,  for  any  given  a  >  1,  we  identify  the  set  K-,  as  the  set  of  all  natural  numbers  k  =  k( a)  such  that 
Sjj:  ^  2 p.  Based  on  the  expression  of  the  terms  Sy<k  given  in  (22),  the  condition  S-,k  ^  2 p,  becomes  equivalent 

to  cos((a+i)— )  ^  _  i  ,  and  furthermore  to  (2j  +  l)n  -  ^  ^  (2k  +  1)  ^  (2j  +  \)n  +  ,  where 

j  =  0, 1, 2, . . . ,  oo.  As  a  result,  we  can  describe  the  set  Ka  as 
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{  k  :  2k  +  1  G 


In 


2n 


(2/  +  1)  7^77  —  1 )  (2y  +  1)  7^77  +  1 


<P{ «) 


(p(a) 


j=  0,1,2,...,  00 


(26) 


For  any  given  a  >  1,  Ky  is  countably  infinite  because  there  is  always  at  least  one  odd  number  of  the  form 
2k  +  1 ,  in  an  interval  of  two-unit  length  such  as  [(2/  +  1)  ^  —  1 ,  (2j  +  1)  ^  +1]  for  j  =  0, 1 , 2, , . . ,  oo. 
This  proves  that  the  terms  Sa^,  and  therefore  the  stress  amplitude,  will  reach  or  exceed  the  value  2 p  for  some 
value  of  k,  for  every  a  >  1.  From  here  it  follows  that  a  design  is  optimal  iff  Sak  =  2 p  for  all  k  e  K( a).  This  is 
equivalent  to  the  fact  that  [2k  +  1  =  (2 j  +  1)  ^  ±  1]  for  j  =  0, 1,2, . . . ,  oo,  which  combined  with  (16),  (17) 
and  (25),  leads  to  the  following  necessary  and  sufficient  condition  for  the  optimal  values  of  the  design 
parameter  a, 


_  1  +  cos  (<?„) 

1  -  COS(<j \)  ’ 


i  ^  n  n 

where  0  <  <pn  =  -  ^  - 
n  2 


and  n  =  2, 3, ,  oo. 


(27) 


Substituting  these  optimal  values  into  the  stress  formulas  Ta.k  and  Sv  k,  we  find  that  these  periodic  functions 
have  a  time-period  of  2 n,  where  the  stress  reaches  its  first  peak  of  2 p  at  the  nth  time  interval.  These  results, 
applied  to  the  time  history  profiles  involving  the  analytical  solutions  (15)  and  (22),  match  with  the  graphical 
output  given  in  Figs.  4a  and  b  and  5a  and  b,  included  in  the  next  page.  The  normalized  values  for  stress  and 
time  used  in  Figs.  4a  and  b  and  5a  and  b,  are  calculated  as,  a*  =  jj,  and  2T  =  (. 

Flere  p,  x  are  the  initially  given  constants  representing  the  stress  loading,  and  the  transit  time  through  the 
elastic  strip.  We  also  notice  that  from  (5),  the  stress  jump  at  the  optimal  designs  with  parameter  a„  takes 
values  [er]T  =  (1  —  cos(jj))[ff]j,  and  [er]R  =  -  cos0[ff]t.  In  general,  for  the  designs  with  parameter, 


1  +  cos  q>j 

<xM  =  t - — ,  where  0  <  q>  -„ 

1  -  cos  q>j  h 


(2y  +  j>  n 


(28) 


there  are  time  intervals  when  the  stress  takes  the  value  2 p.  Here  j  and  n  are  natural  numbers,  chosen  to 
satisfy  the  bounds  for  cpjn  given  in  (28). 

The  following  examples,  provide  a  means  to  better  understand  the  structure  of  the  Ky  set  and  the 
conclusions  made  in  this  section. 


Fig.  4.  Optimal  stress  history  when  x  =  1  and  a.  =  3:  (a)  layer  1:  {  =  0.25.  (b)  layer  2:  {  =  0.75. 
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Fig.  5.  Optimal  stress  history  when  t  =  1  and  a  =  32.163:  (a)  layer  1:  £  =  0.25,  (b)  layer  2:  £  =  0.75. 


Example  1.  a  =  ]^os(ft/3)  =  3-  According  to  (27),  such  a  design  obtained  for  n  =  3,  is  optimal.  In  this  case 
cp(  a)  =|  and  (26)  implies  that, 

K«  =  {k:  2A-  +  1  G  [6(2/  +  1)  -  1, 6(2/  +  1)  +  1],  j  =  0, 1, 2, . . . , oo} 

=  {67  +  2,67  +  3,  j  =  0,1, 2,..., 00}.  (29) 


Based  on  the  structure  of  Ka,  the  stress  amplitude  reaches  its  peaks  when  SXt6j+2  =  2 p,  and  67+3  =  2 p.  As 
predicted  for  n  =  3,  the  first  peak  is  reached  at  the  3rd  time  interval  with  6  units  time-period.  This  can  be 
verified  by  the  graphical  output  given  in  Fig.  4b. 


Example  2.  a  =  Here  tp( a)  =  and  from  (26)  we  obtain, 

7  =  0, 1,2,..., 00  !>. 


»s(57t/  11) 

Ka  =  <j  k  :  2k  +  1  G 


22(27+  1)  |  22(27+  1)  |  { 


(30) 


Although  the  design  parameter  a  is  not  optimal,  it  satisfies  the  condition  (28)  for  j  =2  and  n  =  11.  One  can 
easily  check  that  although  S2  >  2 p  and  S6  >  2 p,  there  are  indices  such  as  {10, 11}  c  Kx  for  which  .S'i0  =  2 p 
and  =  2p,  as  predicted. 

Example  3.  a  =  In  this  case  we  have  that  <p(a)  =  n/s/6  and  from  (26)  we  obtain, 


Ka  =  {k:2k+l  G  [(4/-  +  2)v/6-  1 ,  (4y  +  2) x/6  +  1],  j  =  0, 1, 2, . . . ,  00}. 


(31) 


In  this  case,  none  of  the  conditions  (27)  or  (28)  is  satisfied.  This  means  that  the  design  is  not  optimal  and 
that  all  the  terms  S7k  for  k  €  K.,,  will  exceed  the  value  2 p  without  ever  reaching  it,  i.e.  we  expect  that 
S^k  >  2 p,  Wk  GKa. 

In  conclusion,  in  this  work,  we  solve  the  problem  of  minimizing  the  stress  amplitude  for  a  two-layered 
elastic  strip,  with  layers  of  equal  wave  travel  time,  and  subjected  to  transient  loading.  Our  results  provide 
explicit  formulas  for  the  stress  propagation  and  optimal  designs  in  one-dimension,  and  apply  for  a  two¬ 
layered  elastic  strip  of  unequal  layer  lengths. 
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3.  Numerical  results  using  DYNA3D/GLO  optimization  software 

Since  exact  analytical  solutions  to  optimal  design  problems  involving  transiently  loaded  multilayered 
structures  are  rare  in  the  literature,  analysts  have  relied  on  hybrid  computational  methods  which  link 
formal  nonlinear  parameter  estimation  algorithms  with  computational  finite  element  methods  for  designing 
structures  subjected  to  transient  loadings.  The  results  of  the  previous  section  provide  useful  benchmarks 
which  can  be  used  to  validate  these  hybrid  optimization  design  tools.  Such  a  design  tool  was  recently 
developed  by  Gazonas  and  Randers-Pehrson  (2001),  by  coupling  the  global-local  optimization  software 
package  known  as  GLO  (Murphy,  1999)  with  the  three-dimensional  transient  structural  dynamics  finite 
element  code  known  as  DYNA3D  (Hallquist  and  Whirley,  1989).  The  local  optimization  algorithm  within 
GLO  is  based  upon  the  nonlinear  optimization  program  NLQPEB  that  was  developed  by  Baker  (1992)  for 
optimization  of  shaped  charge  jets. 

The  optimization  software  was  tested  to  see  if  it  could  find  any  of  the  optimal  points  given  by  the  op¬ 
timality  condition  (27).  The  boundary  value  problem  illustrated  in  Fig.  2  was  modeled  using  DYNA3D 
with  60  hexahedral  finite  elements  through  the  thickness.  The  optimization  strategy  is  based  upon  a  global- 
local  approach  since  the  optimal  points  in  layer  2  contain  an  infinite  number  of  minima,  four  of  which  are 
illustrated  in  Fig.  6.  The  plots  in  Fig.  6  were  determined  using  our  analytical  expressions  for  stress,  S^k  and 
TV)k,  through  parametric  variation  of  a  over  the  range  0-10,  and  k  to  a  time  range  out  to  50  units. 

The  global  optimization  uses  a  discrete  scheme  which  simply  subdivides  the  design  parameter  a  into 
uniformly  spaced  regions  from  which  the  local  scheme  can  begin  searching  for  extrema.  The  local  scheme 
seeks  to  minimize  the  maximum  stress  in  layer  2  using  a  BFGS  variable-metric,  sequential-quadric-pro¬ 
gramming  algorithm  (Dennis  and  Schnabel,  1983),  with  a  modified  Powell  merit  function  (Powell,  1977). 
Since  the  number  of  optimal  designs  is  infinite  in  nature,  the  search  was  limited,  rather  arbitrarily,  to  a 
region  involving  the  first  eight  optimal  design  points.  Table  1  shows  that  the  hybrid  DYNA3D/GLO 
optimization  software  successfully  found  all  of  the  first  eight  optimal  design  points  to  a  reasonable  degree 
of  accuracy.  This  result  both  validates  the  DYNA3D/GLO  optimization  software  and  corroborates  the 
optimality  condition  obtained  earlier. 

An  example  of  optimal  layer  design  for  a  two-layered  elastic  strip  of  equal  lengths,  is  a  design  made  of 
tungsten  and  lead  (a  =  3)  or  a  design  made  of  tungsten  and  nylon  (a  =  32.163).  The  proportionality  between 
the  density  and  elastic  modulus  in  both  layers,  provide  a  constant  wave  speed  and  therefore  equal  travel 


a 


Fig.  6.  Maximum  stress  vs.  a  in  a  two-layered  finite  strip  showing  the  first  four  optimal  a  points  given  by  optimality  condition  (27). 
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Table  1 

Comparison  of  analytical  and  DYNA3D/GLO  derived  a  values 


6427 


n 

2 

3 

4 

5 

6 

7 

8 

9 

«  (Eq.  (27)) 

1 

3 

5.828 

9.472 

13.983 

19.196 

25.270 

32.163 

a  (DYNA/GLO) 

1 

3 

5.829 

9.470 

13.924 

19.182 

25.247 

32.120 

times.  These  examples  appear  in  Figs.  4  and  5,  where  it  can  be  seen  that  the  stresses  never  exceed  the  value  of 
2 p,  a  fact  that  qualifies  them  to  be  optimal  designs.  One  can  also  observe  that  the  time  of  onset  of  the 
maximum  stress  in  each  layer  increases  with  the  design  parameter  a.  This  means  that  if  a  design  that  delays 
the  time  of  occurrence  of  the  maximum  stress  is  sought,  then  the  value  of  the  optimal  design  parameter  a„ 
given  in  (27)  should  be  maximized  (within  practical  bounds).  Because  of  such  good  agreement  between  the 
analytical  results  obtained  previously  and  DYNA3D  solutions,  these  findings  support  the  use  of  purely 
computational  means  to  establish  optimal  designs  for  transient  problems  of  this  nature. 

Several  other  problems  with  different  boundary  conditions  were  then  investigated  using  the  DYNA3D/ 
GLO  software,  including  a  two-layered  free-free  system  with  Heaviside  loading,  and  two-layered  free-free 
and  free-fixed  systems  with  rectangular  pulse  loading  p(t)  =  p(H{t)  —H{t—  t0)).  Interestingly,  none  of 
these  boundary-value  problems  provided  optimal  solutions  with  behavior  similar  to  the  optimality  con¬ 
dition  (27).  However,  for  the  free-fixed  problem  subjected  to  the  rectangular  pulse  loading,  as  the  value  of  t0 
increased  and  the  loading  approached  that  of  a  Heaviside  loading,  the  optimal  points  (27)  were  gradually 
recovered. 


4.  Summary 

In  this  paper,  we  investigated  stress  wave  propagation  in  a  two-layered  elastic  strip  subjected  to  free- 
fixed  boundary  conditions  under  a  Heaviside  stress  loading.  Explicit  expressions  for  the  stress  in  each 
layer  were  obtained,  and  the  countably  infinite  set  of  optimal  designs  which  provide  the  smallest  stress 
amplitude  was  identified.  A  time  delay  benefit  when  using  optimal  designs  other  than  the  homogeneous, 
was  discovered.  The  analytical  expressions  for  the  stress  in  each  layer  compared  well  with  the 
stress  history  profile  predicted  when  using  computational  finite  element  methods.  In  addition,  the  op¬ 
timality  condition  (27)  provided  means  to  validate  a  method  based  upon  DYNA3D/GLO  optimization 
software. 

The  validation  of  such  computational  optimization  algorithms  against  known  exact  optimal  solutions 
enables  us  to  consider  a  broader  class  of  impact  and  penetration  problems  using  more  realistic  constitutive 
behaviors  applicable  to  multilayered  armor. 
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